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Abstract 

Electric  power  was  estimated  in  case  of  the  large-scale  flat  thermoelectric  panels  exposed  to  two  thermal  fluids.  The  output  powers 
of  the  proposed  15  systems  were  analytically  deduced  from  heat  transfer  theory,  and  written  by  non-dimensional  functions  to  reflect  the 
characteristics  of  system  design.  The  maximum  output  was  the  largest  in  the  ideal  isothermal  systems.  In  the  other  realistic  systems,  it 
was  the  largest  for  the  system  of  the  counter  flow  with  one  panel.  The  multiplication  of  thermoelectric  panels  can  shorten  significantly  the 
device  area,  although  the  output  from  the  multi-panels  decreases  a  few  percent.  The  worse  heat  conductivity,  X  of  thermoelectric  materials 
and  the  better  heat  transfer  at  the  surfaces  are  desired  for  these  fluid  systems,  in  addition  to  the  better  figure  of  merit,  Z. 
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1.  Introduction 

Seebeck  effect  on  the  junctions  of  two  different  materi¬ 
als  can  generate  the  thermoelectric  power,  P  depending  on 
a  temperature  difference.  Because  we  can  not  use  the  infini- 
tively  large  heat  bath  practically,  this  ideal  case  should  be 
modified.  It  is  designed  that  a  thermal  energy  is  extracted 
by  a  fluid  and  passed  to  the  thermoelectric  generator.  There, 
hot  and  cold  fluids  offer  the  temperature  difference  between 
the  two  surfaces  of  the  plane  thermoelectric  modules.  The 
thermoelectric  motive  force,  E  in  the  open  circuit  is  the  sum 
of  multiplication  of  the  relative  Seebeck  coefficient,  a  and 
the  temperature  difference,  AT  over  all  the  serial  connec¬ 
tions.  The  problem  to  obtain  the  larger  P  is,  therefore,  how 
we  give  the  lager  AT  for  all  the  thermo-elements,  because 
a  is  determined  only  by  the  materials  used  for  the  thermo¬ 
electric  panel. 

The  heat  applied  at  the  hot  junction  is  scattered  and  lost  at 
the  cold  junction,  corresponding  to  the  temperature  gradient 
inside  the  thermoelectric  panel.  The  heat  from  the  hot  fluid 
is  passed  to  the  panel  surface  by  heat  transfer,  to  the  cold  sur¬ 
face  by  heat  conduction  through  the  solid,  and  finally  to  the 
cold  fluid  again  by  heat  transfer.  The  fluids  are  resultantly 
warmed  or  cooled  along  the  flow  path,  and  their  tempera¬ 
ture  profile  through  the  path,  T(x),  changes  as  a  function  of 
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position,  x.  This  kind  of  problem  is  known  as  the  problem  at 
heat  exchanger  through  an  isolator  [1],  Additionally  as  the 
special  feature  in  the  thermoelectric  power  generation,  the 
heat  transfer  due  to  Peltier  effect  and  Joule  heating  should 
be  considered  [2,3]. 

The  purpose  of  this  work  is  to  show  the  mathematical  ex¬ 
pression  of  the  maximum  power  extractible  from  the  ther¬ 
moelectric  panel  which  is  heated  by  the  hot  fluid  and  cooled 
by  the  cold  fluid.  The  heat  transfer  through  the  panel  and  the 
temperature  change  along  the  fluid  path  are  analyzed  math¬ 
ematically.  Mathiprakasam  et  al.  reported  the  precise  anal¬ 
ysis  using  Laplace  transformation  [3].  Their  solution  was 
very  complex  form  due  to  Peltier  effect  and  Joule  heating 
although  a  sheet  of  panel  was  studied. 

It  was  expected  that  a  compact  system  for  the  large-scale 
power  generation  could  be  realized  by  the  multi-sheets  of 
panel,  which  were  efficiently  given  the  larger  AT  by  the 
multi-fluid  paths  [4,5].  As  we  proposed  previously  the  con¬ 
cept  to  use  this  kind  of  multi-layered  panels  [4],  the  heat 
loss  through  a  thermoelectric  panel  is  recovered  by  another 
fluid,  and  the  captured  heat  is  passed  to  another  thermoelec¬ 
tric  panel.  We  hope  that  the  efficient  recovery  of  thermal 
energy  in  a  restricted  space  may  become  possible  by  solv¬ 
ing  the  necessary  thermal  and  fluid  conditions.  Esarte  et  al. 
considered  the  temperature  drop  during  circulation  of  fluid, 
and  showed  that  it  was  fairly  large  [5]. 

Here  we  propose  15  kinds  of  thermoelectric  power 
generation  system  with  multi-panels.  Applying  the  heat 
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transfer  analysis,  we  solve  analytically  their  temperature 
profiles  T(x)  and  output  power,  P  using  the  matrix  analysis. 
Non-dimensional  descriptions  enable  us  to  select  the  best 
system,  independently  of  the  thermophysical  properties  of 
the  thermoelectric  materials  and  working  fluids. 


2.  Basic  assumptions  and  models 

2.1.  Flat  plane  panels 

It  is  assumed  that  all  of  our  thermoelectric  panels  are 
flat,  and  that  they  consist  of  the  thermoelements  with  a 
single  layer  of  11 -type  pn  junctions,  as  illustrated  in  Fig.  1. 
The  thermoelements  are  homogeneously  aligned  perpen¬ 
dicular  to  heat  flow,  combined  tightly  without  open  space, 
and  connected  electrically  in  series.  The  hot  and  cold  fluids 
flow  along  the  both  flat  surfaces  of  the  panel,  and  they  are 
isolated  by  the  panel.  Practically,  as  shown  in  Fig.  1,  there 
exist  the  electrodes  connecting  the  thermoelements,  the 
open  space  between  them,  the  protective  film  between  the 
fluid  and  the  thermoelements,  and  the  fins  that  enhance  the 
heat  exchange.  Esarte  et  al.  proposed  to  circulate  the  fluid 
in  a  path  to  enhance  heat  transfer  [5].  Here  for  simplicity 
we  neglected  these  effects  on  heat  transfer,  or  considered 
them  as  the  heat  transfer  coefficient  at  the  surfaces. 

This  work  will  discuss  the  shape  and  system  holding 
multi-planes.  The  number  of  the  planes,  n  the  directions  of 
hot  and  cold  fluids,  and  the  path  shape  for  fluids  are  classi¬ 
fied  systematically  as  below. 

2.2.  Directions  of  fluids 

Firstly  the  fluid  directions  are  classified  into  four  types, 
depending  on  the  direction  of  cold  fluid  at  the  exist  of  the 
corresponding  hot  fluid:  parallel  flow  type  (P),  counter  flow 
type  (C),  isothermsal  type  (I)  and  the  mixed  type  (X),  as 
shown  in  Fig.  2.  We  assume  only  one  path  for  each  fluid. 
Namely,  we  set  the  boundary  conditions  that  each  system 
has  only  two  paths  connecting  in  series. 

The  parallel  flow  type  is  the  system  that  the  two  fluids 
flow  in  the  same  direction,  between  which  the  thermoelectric 
panel  is  inserted.  The  fluids  flow  in  the  opposite  direction 


Fig.  1.  Simplification  of  thermoelectric  module.  Two  fluids  flow  over  the 
smooth  flat  panel. 


Fig.  2.  Classification  by  the  direction  of  fluid  flow. 


at  the  counter  flow  type.  These  two  types  are  often  selected 
for  heat  exchangers  using  fluids  [1]. 

The  isothermal  type  is  the  extreme  case  that  whole  a 
surface  of  the  most  outside  of  the  system  can  be  kept  at 
a  constant  temperature  by  the  infinitively  large  heat  bath. 
This  models,  for  example,  the  heat  utilization  of  hot  furnace 
wall,  of  phase  transformation  such  as  boiling  [6],  or  of  huge 
amount  of  seawater  as  the  coolant.  Here  we  discuss  the  case 
that  the  cold  surface  is  kept  isothermally,  as  shown  in  Fig.  2. 

The  mixed  type  uses  the  recycled  fluids  at  the  second 
stage,  after  the  hot  and  cold  fluid  is  cooled  and  heated,  re¬ 
spectively,  at  the  first  heat  exchange.  The  increment  of  output 
power  is  expected  by  getting  compulsorily  double  amount 
of  fluid  and  a  medium  temperature  between  those  of  hot  and 
cold  fluids. 

2.5.  Path  shape 

Secondly  the  path  shapes  of  fluids  are  classified  into  three 
types,  depending  on  setting  of  the  boundary  conditions:  me¬ 
andering  type  (M),  helical  type  (H)  and  branched  type  (B), 
as  shown  in  Fig.  3. 

The  meandering  type  bends  the  hot  fluid  paths,  recovers 
the  exhaust  heat  from  the  panel  to  which  the  fluid  gave  the 
heat,  and  will  pass  the  heat  to  the  next  panel.  This  type  aims 
at  the  heat  recycling  and  the  miniaturization  of  the  system. 
Here  we  assume  only  one  path  for  the  cold  fluid. 

The  helical  type  aims  at  the  homogeneous  AT  for  all  the 
panels,  by  flowing  the  two  fluids  in  counter  flow.  The  path 
connection  is  illustrated  as  broken  lines  in  Fig.  3,  however, 
the  seamless  three-dimensional  connection  can  be  realized 
in  the  roll  pancake  model  as  shown  in  Fig.  4(a),  and  in  the 
double  helical  model  as  shown  in  Fig.  4(b).  Their  detailed 
analysis  in  the  cylindrical  coordinate  will  be  reported  sepa¬ 
rately. 

The  branched  type  is  the  system  that  the  two  fluids  sepa¬ 
rate  equally  into  a  few  paths  so  that  the  two  fluids  flow  side 
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Fig.  3.  Classification  by  the  shape  of  fluid  path. 


Fig.  4.  Example  of  three-dimensional  helical  system:  (a)  roll  pancake 
model;  (b)  double  helical  model. 


by  side.  The  number  of  cold  paths  is  taken  larger  than  that  of 
hot  paths,  in  case  of  H  and  B  systems,  because  the  number 
of  paths  becomes  odd  when  the  number  of  panels  is  even. 

Based  on  the  above  classifications,  the  systems  are  named 
after  the  number  of  panel  sheets  (n),  the  direction  of  fluid 
flow  (P,  C,  I  or  X)  and  the  path  shape  (M,  H,  B).  Table  1 
lists  the  names  of  systems  analyzed  here.  Under  this  clas¬ 


sification,  the  previous  conceptual  proposal  [4]  corresponds 
to  the  mixture  of  3CM  and  2X. 


3.  Physical  properties  and  conditions 

The  thermophysical  properties  of  the  thermoelectric  ele¬ 
ments  in  the  literature  scattered  due  to  the  impurities  and 
the  difference  in  preparation  [2],  Table  2  shows  a  typical  set 
of  properties  for  an  identical  sample  [7-10]  for  this  work. 
Tables  3  and  4  show  the  fluid  properties  [11]  and  the  param¬ 
eters  for  thermoelectric  panels.  The  parameters  are  approxi¬ 
mated  as  the  values  at  the  average  temperature  of  7],  and  Tc. 

The  BiTe  semiconductor  is  now  available  in  the  market. 
However,  the  large-scale  power  generation  needs  a  large 
amount  of  materials.  Although  the  thermoelectric  properties 
of  iron-base  alloys  are  not  suitable  still  now,  they  are  attrac¬ 
tive  because  of  the  natural  abundance  and  the  good  electrical 
conductivity  [12], 

The  liquid  is  superior  to  the  gaseous  fluid  as  the  thermal 
carrier  because  of  larger  heat  capacity,  however,  the  gas  has 
the  advantage  in  corrosion  of  the  panels  and  it  can  be  often 
used  as  the  exhaust  gases  from  the  industries.  This  study 
assumes  the  usage  of  inert  N2  gas  as  heat  carrier. 

In  general,  heat  transfer  coefficient,  h  depends  mainly  on 
the  hydrodynamic  behavior  of  the  fluid  and  the  shape  of  the 


Table  1 

Classification  and  abbreviation  of  studied  system 


Shape  of  fluid  path 

Direction  of  fluid  flow 

Parallel  (P) 

Counter  (C) 

Mixed  (X) 

Isothermal  (I) 

One  panel 

- 

IP 

1C 

- 

11 

Meandering  (M) 

2PM 

2CM 

- 

21M 

3CM 

4CM 

Multi-panel 

Helical  (H) 

- 

2CH 

- 

- 

3CH 

4CH 

Branched  (B) 

- 

3CB 

- 

- 

4CB 

Others 

- 

- 

- 

2X 

21 
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Table  2 


Specific  values  of  thermoelectric  elements  used  here 


Materials 

Seebeck  coeficient, 
of  (jxV/K) 

Resistivity,  p 
(|ji£2  m) 

Thermal  conductivity, 

X  (W/Km) 

Figure  of  merit,  Z 
(l/K) 

Fe  23  at.%  A1  20  at.%  Si  (p) 

40.35 

1.923 

12.44 

4.483  x  10”5 

Fe  22  at%  A1  in) 

-21.56 

1.176 

16.13 

Bi  54.3  at.%  Te  (p) 

162 

5.5 

2.06 

2.605  x  10"3 

Bi  64.5  at.%  Te  (n) 

-240 

10 

2.02 

Table  3 

Parameters  for  thermoelectric  power  generation  system 


Variables 

Values  used  in  this  work 

Thermoelectric  device 

Width,  w 

lm 

Length,  / 

Variable 

Thickness,  d 

0.05  m 

Thermal  sources 

Hot  source 

N2  gas  (inlet  Tf  =  500  K) 

Cold  source 

N2  gas  (inlet  Tf  =  300  K) 

Gas  properties 

Specific  heat,  Cp 

1044J/kgK  (at  400  K,  0.1  Pa) 

Mass  flow  rate,  M 

1.00  kg/s 

Heat  transfer  coefficient,  h 

500  W/m2  K 

tube  surface.  It  was  evaluated  as  28.3  W/m2  K,  assuming  that 
the  turbulent  flow,  which  obeys  with  Kays  equation  [13], 
flows  under  the  hydrodynamic  conditions  listed  in  Table  4. 
The  fins  on  the  surfaces  can  improve  the  heat  transfer  coef¬ 
ficient,  h,  to  a  few  to  ten  times  larger  values  [14].  This  study 
sets  h  as  500  W/m2  K. 


4.  Equations  for  output  power 


Fig.  5.  Module  with  n  panels.  Direction  of  heat  flow  among  fluids  is 
illustrated. 


equations  can  be  written  from  the  heat  transfers  through 
panels: 


/VTiCp.i 


d71(x) 

dx 


=  TKiw(Jx{x)-T2{x)) 


0) 


4.1.  Deduction  of  equations 

Considering  the  system  consisted  of  n  sheets  of  panels 
as  shown  in  Fig.  5,  a  set  of  n  +  1  simultaneous  derivative 


d7/(.r) 

MjCpj — j—  =  ±Ki-lw(Ti-l{x)  -  Ti(x )) 
dx 

=FK,  w(7Kx)  -  Ti+  i(x))  (1  <  i  <  n) 

(2) 


Table  4 

Hydrodynamic  and  boundary  conditions 

System  Hydrodynamic  conditions  Boundary  conditions 


IP 

Mi 

> 

0,  m2 

>  0 

1C 

Mi 

> 

0,  m2 

<  0 

11 

Mi 

> 

0 

2PM 

Mi 

> 

0,  m2 

=  — M  i ,  A^3 

<  0 

2CM 

Mi 

> 

0,  m2 

=  — M[ ,  Mf, 

>  0 

2CH 

Mi 

> 

0,  m2 

A 

o 

II 

Mi 

2IM 

Mi 

> 

0,  m2 

=  -Ml 

21 

Mi 

> 

0,  m2 

=  -Ml 

2X 

Mi 

> 

0,  m2 

=  -(Ml  +  Mi),  Mi  >  0 

3  CM 

Mi 

> 

0,  m2 

=  -Mi .  Mi 

=  Mi ,  M4  >0 

3CH 

Mi 

> 

0,  m2 

=  M4,  Mi  = 

=  Mi ,  M4  <0 

3CB 

Mi 

> 

0,  m2 

=  M4,  Mi  = 

=  Mi ,  M4  <  0 

4CM 

Mi 

> 

0,  m2 

=  — M  i ,  Mf, 

=  Mi ,  M4  =  —Mj 

,  M5  >  0 

4CH 

Mi 

V 

p 

J; 

=  M4,  M3  = 

M\,  Mt,  <  0,  M5  : 

=  Mi 

4CB 

Mi 

> 

0,  m2 

=  m4,  m3  = 

-  Mi,  M4  <  0.  M5 

=  Mi 

7i(0)  =  Tf,  T2( 0)  =  Tf 

T,(0)  =  7^",  T2(t)  =  Tf 

THO)  =  Tf,  T2(x)  =  If  for  all  x 

7i(0)  =  Tf,  Ti(l)  =  TiC/),  73®  =  Tf 

7-1(0)  =  Tf,  ri©  =  72(0,  73(0)  =  Tf 

7-1  ( 0)  =  Tf,  7-1(0  =  7-3(0),  7-2(0  =  Tf 

7-i(0)  =  Tf,  TOO  =  7-2(0,  7-3  (x)  =  Tf  for  all  * 

T-i  (0)  =  Tf,  7-2  (0  =  Tf,  T3(x)  =  Tf  for  all  x 

7-1(0)  =  Tf,  72(0  =  (Mi 7-i(0  +  M37’3(0)/(M1  +  M3),  71,(0)  =  Tf 

7-i(0)  =  Tf,  7-1(0  =  72(0  ,  72(0)  =  73(0),  7i(0  =  Tf 

7-i(0)  =  Tf,  7-i(0  =  73(0),  72©  =  74(0),  r4(0  =  Tf 

7-i(0)  =  Tf,  7-2 (0  =  Tf,  7-3(0)  =  Tf ,  T4(l)  =  Tf 

7-1(0)  =  Tf,  7i(0  =  72(0  ,  7-2(0)  =  73(0),  T3(l)  =  T4(l),  Ts(0)  =  Tf 

7-1(0)  =  Tf,  7-1©  =  73(0),  72©  =  74(0),  73(0  =  75(0),  74(0  =  7f 

7-i(0)  =  Tf,  7-i(0  =  Tf ,  7-2 (0)  =  Tf,  73(0  =  7f ,  74(0)  =  T;1" 
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dP,I+i(x) 

M„+1Cp,„+i  "+1  =  ±Knw(Tn(x)  -  r„+1(jt))  (3) 

dx 

where  M,-,  Cp.(,  7/  are  mass  flow  rate,  heat  capacity  and  tem¬ 
perature  for  fluid  i,  w  the  width  of  the  panel.  The  homoge¬ 
neous  temperature  profile  along  y  axis  is  assumed  inside  the 
turbulent  flow.  The  choice  of  ±  depends  on  the  path  con¬ 
dition,  Kj  the  over-all  heat  transfer  coefficient  through  the 
panel  i, 

Ki  =  - - -  (4) 

(1/Ah./)  +  (d/k)  +  (l//»c.i) 


hhj  and  hcj  are  the  heat  transfer  coefficients  between  hot 
fluid  and  the  panel  and  between  cold  fluid  and  the  panel, 
respectively,  d  and  X  the  thickness  of  the  panel  and  the 
average  heat  conductivity  in  the  panel,  respectively. 

.  _  7.pSp  "F 

k  =  - ; -  D) 

Sp  +  sn 


where  a  is  the  difference  between  those  of  p  and  n  elements, 
nx  and  n  v  the  number  densities  of  thermoelectric  pairs  in  the 
direction  x  and  y,  respectively.  Note  that  nxny  =  1  /(.vp  +  ,vn). 
The  electric  resistance,  P,  is  given  as: 

Ri  =  nxtiy  wld  (  —  +  —  )  (12) 

V  sp  ■Sn  / 


where  pp  and  pn  are  the  specific  resistivity  of  p  and  n 
elements.  The  maximum  output  power,  P  optimized  by 
balancing  the  internal  and  external  resistance  as  used  in 
normal  way: 

E2 


P  = 


4  £?=,*/ 


(13) 


This  work  uses  hereafter  this  P,  where  the  internal  and 
external  resistance  are  the  same. 


4.2.  Output  power  of  IP  system 


where  A.p  and  are  heat  conductivity  of  p  and  n  type  el¬ 
ements,  respectively,  sp  and  sn  the  cross-sectional  area  of 
these  elements.  For  simplicity,  heat  transfer  by  Peltier  ef¬ 
fect  and  exothermal  heat  by  Joule  effect  are  neglected.  The 
detailed  considerations  on  this  simplicity  will  be  reported 
separately. 

By  replacing: 


K  = 


KjW 

MjCpj 


(6) 


the  simultaneous  equations  can  be  written  in  the  matrix  form 


The  simultaneous  equations  are  solved  at  the  IP  system, 
where  one  sheet  of  panel  is  exposed  in  parallel  flow. 


7hU)  = 


Tc(x)  = 


MhCp,hr“  +  McCp,c((P‘n  -  Tf  ) 
xexp  i  —DKwx)  +  ) 

H^hF'p.h  +  HfcCp,  c 

Pcm)  exp  (—DKwx)) 


MhCp,h(K  ~ 

+McCpxT'c 


(T 


A^hCp.h  +  McCj 


P.C 


(14) 


(15) 


as, 

'  T\  (x) 

”  =FAn 

±An 

0 

'  T\  (x) 

d 

Ti(x) 

±An 

T(Al2  +  A22) 

±^22 

T2(x) 

dx 

Tn  (x) 

i-^72— l,n  — 1,«  “1“  An,n) 

i-^72,72 

Tn(x) 

_  Tn+ i(x)  _ 

_0 

72, 72  +  1 

T  ^72, 72+1 

_  Tn+ i(x)  _ 

The  equation  is  solved  for  x  by  setting  the  boundary  condi¬ 
tions,  as  listed  in  Table  4.  The  mathematical  software,  Math- 
ematica,  was  partially  used  as  a  tool  [15].  The  solutions, 
7',(x),  are  described  later.  The  temperature  Ojj(x)  at  the  sur¬ 
face  of  the  panel  i  facing  to  the  fluid  j  is  written  using  Tf  x). 

dU(x)  =  Tfx)  -  fi(7)(x)  -  Ti+i(x))  (8) 

m 

eUi+ 1  (x)  =  Ti+1  (x)  +  t^(7)(x)  -  Ti+l  (X))  (9) 

«i+i 

The  electromotive  force,  E,  is  the  summation  of  the  temper¬ 
ature  difference  over  all  the  panels. 


where 


D 


Mh  f  'p.  h 


1 

H7cCpic 


(16) 


The  electromotive  force  is  calculated  by  integrating  the 
above  results  in  Eq.  (11), 


MhMcCpMCp,cnxny&a(T™-T™)^  ~  exp(-DKwl)) 

(MhCp,h  +  McCp,c)X 


(17) 


The  output  power  Pi  is  obtained  using  Eq.  (13). 


n 


£  =  E 

y KiQu  -  &u+ 1) 

1=1 

panel 

II 

rw  nl 

nxny  /  /  a(0i  i(x)  -  0u+]  (x))  dxdy 

0  JO 


(10) 

(ii) 


MlM2C2p  hclcnxnyda2(4-  -  T?)2 

p  _  x{l  -  exp (-DKwl)}2 _ 

1  4wl(MhCp,h  +  McCp,c)2(PP/sP  +  PnAnU2 

Pi  agreed  with  the  solution  reduced  by  Laplace  transfor¬ 
mation  [3],  when  our  assumptions  were  introduced.  The 
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Fig.  6.  Non-dimensional  output  of  one-panel  systems. 


ference.  ip(A)max  is  constant  independent  of  the  materi¬ 
als  and  fluid  properties,  and  it  is  also  the  parameter  deter¬ 
mined  by  the  system  design.  Therefore,  we  define  y(lP)  = 
<£ip(X)max  as  the  system  parameter,  and  <5(1P)  =  0.6282 
as  the  characteristic  length,  specialized  for  IP  system.  These 
are  listed  in  Table  6. 

The  contributions  of  sp  and  ,vn  through  the  terms  of  K  and 
X  were  neglected  because  they  were  small. 

4.3.  Output  power  of  1C  and  II  system 

The  simultaneous  Eq.  (7)  are  solved  at  the  1C  system, 
and  the  thermoelectric  force  E  and  the  output  power  Pi  are 
calculated  in  the  same  procedure  as  described  above. 


temperature  difference,  0h(x)  ~  &c(x),  agreed  with  the  pre¬ 
vious  analysis  of  parallel  flow  [5]. 

Especially  when  all  the  terms  of  M;CPi,  are  equivalent 
and  all  Kj  are  equal  to  a  constant  K .  we  obtain  the  output 
power  P2  in  a  simple  form, 

Kd 

Pi  =  MCp—Z(T™  -  7f  )2L(sp,  sn)0 ip(X)  (19) 

A 


where  Z  is  the  figure  of  merit,  known  as  the  material  param¬ 
eter  [2,10],  and  <J>ip(A)  the  new  non-dimensional  function. 

(ap  -  oin)2 


Z  = 


(•y/ Pp^-p  T"  s/ Plp-n  P 


(20) 


MhMcCp,hCp,c«t»vda(Pi;n  -  7^X1  -  exp (-DKwl)) 
(MhCp^exp(-DKwl)  +  McCp,c)X 

(27) 

MZM?Cl,An*nyda2(T"  -  Tf  )2 

x{l  ~  exp  (-DKwl)}2 _ 

4wl{MhCp,hexp(-DKwl)  +  McCp,c}2 

x(Pp/sp  +  Pn/Sn)X2 


When  MhCp,h  =  |McCp,c|, 

Kd 

P2  =  MCp  —  Z(T™  -  Tf'fLisp,  sn)<Plc{X)  (29) 

A 


Pphp  + 

L(Sp,  5n)  =  - 

(Ap^p  +  A.niS'n)(Pp5'n  +  Pn^p) 

(21) 

(l-e“2X)2 
<PtP(*)  =  lfiz 

(22) 

where  L(sp,  sn)  is  identical  with  Eq.  (21).  The  function  of 
&ic(X)  is  given  in  Table  5,  and  the  characteristic  parameters, 
<5(1C)  =  1  and  <Pic(W)max  =  y(lC)  =  1/16,  are  listed  in 
Table  6. 

Similarly  Pi  for  II  system  is  given  by: 


where  the  non-dimensional  length  of  the  system  X  is  defined 
as. 


Kwl 

MCp 


(23) 


(p\\>(X)  is  characteristic  for  each  system  as  shown  in  Fig.  6 
and  Table  5.  L(sp,  ,vn)  the  function  of  the  cross-sectional 
area,  and  its  maximum  is  1  when: 


Stp 


(24) 


The  maximum  <Pip(X)  =  5.091  x  10  2  when: 

4X  +  1  =  e2X  (25) 


This  Eq.  (25)  can  be  numerically  solved  as  X  =  0.6282. 
P3  is  obtained  after  these  maximizations  of  L(sp,  ,vn)  and 
<Pip(X)  for  Eq.  (19). 

P3  =  <P[?(X)m.diXMCp1^-Z(T™  -  Pcin)2  (26) 

A 

P3  depends  on  the  mass  flow  through  M,  h\t  and  hc,  as  well 
as  the  figure  of  merit  and  the  square  of  temperature  dif- 


Pi  = 


MlClMnxnyaa2{T™-T™)2 
4wl(pp/sp  +  Pn/sn)X2 


1  —  exp 


Kwl 
47hCp  h 


(30) 


It  is  mathematically  proved  that  Pi  for  II  system  is  a  limit 
of  those  for  IP  and  1C  systems,  as  Mc  approaches  infinity. 
P2  for  II  system  is  also  written  in  the  similar  equation  with 
Eq.  (29).  0\i(X),  5(11)  and  y(  1 C )  are  shown  in  Tables  5  and 
6.  As  Mh  approaches  infinity  in  Eq.  (30),  the  system  corre¬ 
sponds  to  the  case  where  the  both  surfaces  are  heated  and 
cooled  by  the  infinitively  large  heat  baths.  We  obtain  the 
well-known  output, 

p  __  nxnywlot2(Th  -  Tc )2 
4d(pp/sp  +  Pn/Sn) 


where  P  becomes  the  maximum  at  i'n/^'p  =  (Pn/Pp)1^2- 
In  a  short  summary,  the  non-dimensional  output  <t>  de¬ 
pends  on  the  fluid  direction,  but  the  formula  can  be  com¬ 
monly  written  only  by  the  same  length  parameter  X.  Three 
0(x)  for  one  sheet  of  panel  are  compared  in  Fig.  6. 
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Table  5 

Calculated  non-dimensional  output  power  0tX) 


System 

IP 

1C 

II 

2PM 

2CM 

2CH 

2IM 

21 

2X 

3CM 

3CH 

3CB 

4CM 

4CH 

4CB 


MX) 

MX) 

MX) 

MX) 

MX) 

MX) 

MX) 

MX) 

MX) 

MX) 

MX) 

MX) 

MX) 

MX) 

MX) 

MX) 


(1  -  e“2x)2 
- s 

16X 

X 

4(1  +  X)2 

(1-e-*)2 

4X 

{  _  i  +  ^2  -  (1  +  J2)^-X  +  2e<1+'/S>x}2 
2{2  —  \fl  +  (2  +  v'2)e2'/2X)2X 
(l-e"x)2 
8X 

1 

2X(2  +  cothX)2 

{— \/5ex/2  +  a/5  cosh  (V5X/2)  +  3  sinh  (V5X/2)}2 
8X{2  +  3  cosh  (a/5X)  +  a/5  sinh  (v^X)} 
e(-i-V5)X{(_3  +  ^)ex,2  _  2V5eM5X/2  +  (3  + 

32XJ-2  +  7  cosh  (a/5X)  +  3a/5  sinh  (v^X)} 

(1  -  e_x)2 
8X 

{- 1  +  (1  +  2X)  cosh  (s/2X))2 
6X{a/2(2  +  X)  cosh  (V2X)  +  sinh  (V2X)}2 

{2X  cosh  (X/y/2)  +  V2  sinh  (X/V^)}2 
3X{(2  +  4X)  cosh  IX/ a/2)  +  -J2( 2  +  X)  sinh  (X/a/2)}2 

[V2{— 1  +  (1  +  4X)  cosh  (2V2X)}  +  (1  +  6X)  sinh  (2s/2X)]2 
12X[V2{-1  +  (3  +  4X)  cosh  (2s/2X)}  +  2(2  +  3X)  sinh  (2V^X)]2 

[— 2ea+3V5W2(— 1  +  ex)(ll  +3ex)  +  ex{-ll  _  5^/5  +  (7  —  3V5)ex)  +  e(1+3v^)X{— 11  -  5a/5  +  (7  +  Sv^e*} 
+e'/3X[10  _  6^/5  +  ex{- 1 1  +  3V^  +  (-3  +  y,5)ex}]  +  e2^!- 10  -  6^  +  ex{l  1  +  3^5  +  (3  +  V,5)ex}] 

1  +e(i+V5)X/2|-_  i  _  ^5  +  ex{8  _  6^5  -  (7  -  3y,5)exi]  +  e(1+5'/3>x/2[-l  +  ^5  +  ex{8  +  6s/5  -  (7  +  3y,5)ex}]]2 

4X  {(5  _  3y/5)ex  +  10e^x  +  (5  +  3A/5)e(1+2^x)2{(-3  +  ex)(-l  +  e^x)  +  V5(-l  +  ex)(l  +  e^x)}2 

{/5e -  5V5e<‘+^>x  -  5e2<1+^>x  +  5-v/5e<1+3^>x  -  2(5  +  V5)e3«+^>x/2  -  2(-5  +  V5)^+3^x'2 
1  +2(5  +  y5)e<‘+5'/W2  +  2( — 5  +  V5)e<3+W5)x/2}2 

4X  [15eVsx  _  8eU+V5)X  _  5e(2+V5)X  _  20e3(i+V5)X/2  _  20e(3W5)X/2  +  4(5  _  2V5)e<1+^x/2 

+(9  -  10^)ex  +  4(5  +  2v/5)e(1+3'/5)X/2  +  (9  +  10v/5)e(l+2^5)x)2 

{  -  1  +e5x(l  +45X))2 
4X{2  —  9e5X(3  +  10X))2 


The  outputs  are  calculated  using  the  materials  parameters 
in  Table  2  for  the  BiTe  semiconductors  and  for  the  FeAl  al¬ 
loys.  As  shown  in  Fig.  7,  the  maximum  power  for  BiTe  semi¬ 
conductor  requires  a  characteristic  panel  length  3.76  times 
longer  than  that  of  FeAl  alloys,  but  it  gains  106  times  larger 
output  power,  when  the  panel  thickness  d  is  the  same.  Be¬ 
cause  the  ratio  of  Z  in  both  materials  is  58.1,  the  worse  heat 
conductance  of  BiTe  causes  the  larger  difference  in  output 
power.  For  the  materials  selection  for  thermoelectric  panels, 
therefore,  the  total  output  power  and  the  prices  of  materials 
and  fluid  will  be  considered  entirely,  as  we  pointed  previ¬ 
ously  [12].  In  case  of  thinner  BiTe  panel  ( d  —  0.005  m), 
for  example,  the  output  power  decreased  although  the  panel 
length  can  be  shortened.  This  shows  that  the  optimal  thick¬ 
ness  of  the  panel  exists  depending  on  a  certain  length,  as  re- 


Fig.  7.  Output  power  of  1C  system,  where  N2  gases  of  1  kg/s  at  500  and 
300  K  are  flown  over  a  panel  of  1  m  width  with  counter  flow,  d  is  the 
panel  thickness. 


208 


R.O.  Suzuki,  D.  Tanaka/ Journal  of  Power  Sources  122  (2003)  201-209 


Table  6 

Analyzed  non-dimensional  parameters 


System 

System 
parameter,  y 

Characteristic 
length,  5 

Output  per  unit 
length,  y/8/n 

ip 

0.050908 

0.62822 

0.08104 

1C 

1/16  =  0.0625 

1 

0.06250 

11 

0.101816 

1.25643 

0.08104 

2PM 

0.050881 

1.77891 

0.01430 

2CM 

0.050908 

1.25643 

0.02026 

2CH 

0.058289 

0.41581 

0.07009 

2X 

0.050908 

1.25643 

0.02026 

21 

0.091863 

0.63631 

0.07218 

2IM 

0.092222 

1 .69764 

0.02716 

3  CM 

0.045203 

1.73604 

0.00858 

3CH 

0.059231 

0.28820 

0.06851 

3CB 

0.056444 

0.30489 

0.06171 

4CM 

0.040966 

2.07989 

0.00492 

4CH 

0.060183 

0.22504 

0.06686 

4CB 

0.059793 

0.24565 

0.06085 

ported  experimentally  [16].  Fig.  8  shows  the  optimal  thick¬ 
ness  evaluated  from  Eq.  (29).  Our  deduced  equations  are 
thus  the  base  for  system  design  of  thermoelectric  power  gen¬ 
eration. 

4.4.  Output  power  of  n  sheets  of  panels 

The  equations  are  calculated  in  the  system  having  n  sheets 
of  panels.  The  analytical  solutions,  P\,  are  obtained  in  the 
same  sequence  of  Section  4. 1  and  under  the  conditions  listed 
in  Table  4.  The  obtained  equations  are  complex  and  very 
long  expression.  When  they  are  simplified  by  setting  all  K, 
identical  and  M^Cph  —  |McCPiC|,  the  output  power.  If.  in 
all  the  systems  can  be  commonly  described  as  Eq.  (19)  or 
(29)  by  the  non-dimensional  output  <P{X).  They  are  listed  in 
Table  5,  where  the  comparison  among  the  systems  becomes 
possible,  independently  of  thermoelectric  materials  and  of 
the  thermophysical  properties  of  fluids. 

Figs.  9  and  10  show  X  dependency  of  Only  one 
maximum  commonly  exists  for  each  The  system 


Fig.  8.  Output  power  of  1C  system,  where  N2  gases  of  1  kg/s  at  500  and 
300  K  are  flown  over  a  BiTe  panel  of  1  m  width  with  counter  flow. 


Fig.  9.  Non-dimensional  output  of  the  two-panel  systems. 


Fig.  10.  Non-dimensional  output  of  the  three-  and  four-panel  systems. 


parameter  y  and  the  characteristic  length  <5  are  evaluated 
numerically  in  Table  6. 


5.  Discussions 

The  system  parameter,  y  are  summarized  in  Fig.  11. 
Firstly,  /(II),  /(2I)  and  /(2IM)  are  significantly  high,  be¬ 
cause  a  cold  surface  is  kept  at  a  constant  temperature  for 
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these  systems.  It  was  basically  assumed  in  the  isothermal 
type  (I)  that  the  infinitive  fast  heat  supply  and  heat  compen¬ 
sation  was  available  to  keep  a  constant  temperature.  This 
assumption  was  mathematically  realized  that  a  limit  of  IP 
and  1C  system  was  II  system  as  M  approached  to  the  in¬ 
finitive.  The  infinitive  fast  cooling  in  I  type  is  the  reason  of 
high  power. 

Secondly,  the  power  of  1C  is  larger  than  that  of  IP,  while 
the  characteristic  length  of  1C  is  longer.  As  shown  in  Fig.  6, 
<J(1C)  is  larger  than  <J>(1P)  for  all  X.  The  tendency  that  the 
counter  flow  can  generate  larger  power  than  the  parallel  flow 
is  valid  also  in  two  sheets  of  panels:  y( 2CM)  >  y( 2PM). 
However,  the  output  per  unit  length  of  panel  is  the  highest 
for  IP  and  II,  as  listed  in  Table  6.  The  efficiency  in  thermo¬ 
electric  panel  is  better  for  IP  and  II  than  for  1C,  although 
the  efficiency  of  heat  is  the  best  for  1C. 

When  we  compare  by  the  path  shapes  (H,  B  and  M),  the 
helical  type  has  the  highest  y  and  the  shortest  length  8.  The 
branched  type  is  a  little  smaller  y  and  a  little  longer  8  than 
those  of  H  type.  The  geometrical  construction  of  H  type 
seems  possible  as  shown  in  Fig.  4,  although  the  total  path 
length  for  fluid  becomes  longer,  which  causes  larger  pressure 
loss.  Considering  the  pumping  power  of  fluid  [5],  therefore, 
the  latter  may  have  still  the  potential  for  application. 

As  the  number  of  sheets  of  panels  increases  for  H  and 
B  types,  their  system  parameters  y  increase  and  their  char¬ 
acteristic  lengths  8  become  shorter.  In  H  and  B  types,  the 
multi-panel  construction  thus  improves  the  thermoelectric 
properties,  however,  it  decreases  them  for  M  types.  Mean¬ 
dering  type  needs  the  longest  8  in  these  three  types.  Namely, 
heat  recycling  by  M  type  needs  plenty  of  thermo-elements.  If 
we  can  use  the  cheaper  thermoelements  in  future,  the  exhaust 
heat  from  the  thermoelectric  generators  will  be  recovered. 

Except  for  isothermal  types,  the  highest  output  parameter 
y  is  obtained  in  1C  system,  y  becomes  gradually  smaller 
at  4CH,  4CB,  3CH,  2CH,  3CB  systems  in  this  order, 
yfmulti-panels  systems)  gets  closer  to  y(  1 C),  as  the  number 
of  sheets  increases  in  H  and  B  types. 

The  shortest  length  is  achieved  at  4CH  system,  where  the 
difference  with  output  of  1C  system  is  only  3.71%  and  the 
length  of  4CH  is  77.5%  shorter.  Therefore,  the  output  per 
unit  length  of  panel  is  higher  for  H  types  than  1C. 

In  summary,  the  multiplication  of  panels  decreases  the 
output  a  little,  but  it  can  shorten  the  length  of  thermoelectric 
power  generator. 

6.  Conclusion 

This  work  deduced  analytically  the  output  Pi  for  15  sys¬ 
tems  of  flat  thermoelectric  panels,  based  on  heat  transfer 
theory.  Setting  the  heat  permeability  K,  mass  flow  rate  M 
and  specific  heat  of  fluids  Cp  identical,  the  output  P3  can  be 
written  by  non-dimensional  functions  <P(X).  The  maximum 
output  P 3  is  determined  as: 

KrJ 

p3  -  '/.V7CP  — -  z<  y;;,:  -  iff 

A 


The  system  parameter  y  is  characteristic  for  geometry  of 
path  and  the  direction  of  fluid,  y  were  significantly  large 
for  isothermal  systems.  In  the  other  systems,  y  was  largest 
for  the  system  of  the  counter  flow  with  one  panel  (1C).  The 
multiplication  of  thermoelectric  panels  can  shorten  signifi¬ 
cantly  the  system  length,  i.e.  the  device  area,  although  the 
outputs  from  the  multi-panels  decrease  a  few  percent. 

The  effect  of  thermoelectric  material  properties  on  the 
output  in  thermoelectric  power  generation  has  been  consid¬ 
ered  representative  only  by  the  figure  of  merit,  Z.  However, 
as  written  in  Eq.  (32),  X  has  the  contribution  to  the  output 
stronger  than  as  expected  previously,  when  the  fluids  are 
working.  Under  the  fixed  conditions  given  by  such  as  M, 
Cp,  Z  and  AT,  the  output  is  a  function  of  the  term  K/X  in 
Eq.  (32).  The  better  heat  transfer  to/from  the  surface  and  the 
worse  heat  conductivity  of  the  panel  are  desired. 
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